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Abstract 

We study the propagation of null rays and massless fields in a black hole fluctu- 
ating geometry. The metric fluctuations are induced by a small oscillating incoming 
flux of energy. The flux also induces black hole mass oscillations around its aver- 
age value. We assume that the metric fluctuations are described by a statistical 
ensemble. The stochastic variables are the phases and the amplitudes of Fourier 
modes of the fluctuations. By averaging over these variables, we obtain an effective 
propagation for massless fields which is characterized by a critical length defined 
by the amplitude of the metric fluctuations: Smooth wave packets with respect to 
this length are not significantly affected when they are propagated forward in time. 
Concomitantly, we find that the asymptotic properties of Hawking radiation are not 
severely modified. However, backward propagated wave packets are dissipated by 
the metric fluctuations once their blue shifted frequency reaches the inverse critical 
length. All these properties bear many resemblences with those obtained in models 
for black hole radiation based on a modified dispersion relation. This strongly sug- 
gests that the physical origin of these models, which were introduced to confront the 
trans-Planckian problem, comes from the fluctuations of the black hole geometry. 
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1 Introduction 



In his original derivation of black hole radiance, Hawking [0 considered the propagation 
of a linear quantized field in a classical background geometry, that of a collapsing body. 
In this framework, one neglects two effects. One first neglects the gravitational back 
reaction effects, i.e. the consequences of the quantum response of the geometry to the 
energy density of the radiation field. To compute these effects is at present out of reach 
as it requires a better understanding of quantum gravity. 

One also neglects the fluctuations of the geometry which are not due to the energy den- 
sity of the radiation field. Besides quantum mechanical fluctuations of the gravitational 
field itself, there exist also metric fluctuations induced by the quantum fluctuations of 
other fields. The latter can be approximatively described by introducing stochastic noice 
sources in the right-hand side of the Einstein equations [0, |^, ^, |]. In this description, 
one therefore deals with a stochastic ensemble of fluctuating geometries. Our aim is to 
study the propagation of a massless field in such an ensemble. 

To describe metric fluctuations near the black hole horizon we shall use a model 
similar to that considered by York [Q. It is based on the hypothesis that the metric 
fluctuations are driven by a small oscillating flux of energy of an infalling null fluid. In 
our model, the metric fluctuations are represented by a linear superposition with different 
frequencies. Stochasticity will come into the picture by assuming that the amplitudes and 
phases of each mode are stochastic variables. Therefore to obtain the expectation value 
of any observable will require averaging over these variables. As we shall still neglect 
gravitational back reaction effects, we shall still have a linear field equation. In a former 
paper ^ , we analyzed the propagation of a test field in a metric characterized by a given 
classical fluctuation. Therefore, the novelty of the present work is to take an ensemble of 
such fluctuations. 

The influence of black hole metric fluctuations on physical effects in the black hole 
geometry is an interesting and open problem. As we shall demonstrate in the paper, 
one of the justifications of considering a statistical ensemble of metric fluctuations is to 
confront the trans-Planckian problem^], ^ |10|, 0, |12|, 0, |15|, [l^, ll^. This problem 



comes from the fact that low energy modes reaching J'^ emerge from field configurations 
which possess, near the event horizon, arbitrarily high frequencies (as measured in a freely 
falling frame). Since gravitational interactions grow with the energy, one must call into 
question the validity of describing the propagation of these field configurations by free 
field theory. In particular, one might wonder if the low energy asymptotic properties of 
Hawking radiation, namely stationarity and thermality with a temperature determined 
by the surface gravity, are sensitive to the high energy behavior of the theory. 

To answer this point, Unruh made an interesting proposal p, |15]. He considered the 
propagation of sonic waves in an acoustic geometry governed by a wave equation whose 
dispersion relation is bended at high frequency. Indeed, only low frequency phonons 
propagate freely with a given velocity. For frequencies higher than a critical frequency Uc, 
the dispersion relation is no longer linear and dissipation may occur. He then showed that 
the modification of the dispersion relation in no way affects the asymptotic properties of 
Hawking radiation as long as Uc is much bigger than the surface gravity k. The reason 



is that there is an adiabatic decoupling between these two energy scales. In |T^, it was 
conjectured that light propagation near a black hole horizon should also be described by 
an effective mutilated theory and an alternative model "which may be more appropriate to 
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black hole physics" was proposed. The basic argument is that when a Hawking quantum 
is propagated backward close to the horizon, it will interact with a reservoir of modes, 
e.g. the high angular momentum modes which do not reach spatial infinity, or the metric 
fiuctuations induced by these modes 



In this paper, we consider a stochastic ensemble[] of metric fiuctuations to describe 
these interactions. We shall show that light propagation in a stochastic metric indeed 
leads to an effective truncated theory near the event horizon. More precisely, we obtain 
the following. First the critical lenght uj~^ is determined by the amplitude of the metric 
fiuctuations. Secondly, as far as forward propagation is concerned, the evolution of smooth 
wave packets (where smooth means that their in-frequency content is much below Uc) is 
affected only slightly by the metric fluctuations. Thirdly, backward in time propagation 
of wave packets representing Hawking quanta is dramatically modifled only when the blue 
shift factor brings their frequency close to Uc- In this regime, the amplitude of the wave 
packet is rapidely dissipated (backward in time!). 

At this point, the reader might wonder about the physical validity of these results 
since we do not know the precise nature of the metric fluctuations near the event horizon. 
In this respect, the following point should be emphasized. We are not studying the 
fluctuations themselves but only their effects on light propagation. As argued by Feynman 



and Hibbs|^, the effective propagation obtained by tracing out the degrees of freedom 
of the environment does not depend on the precise nature of the interactions of the 
test particle with it. For this universality to apply, one should neither ask too detailed 
questions (e.g. during too small time intervals) nor consider too strong interactions leading 
to signiflcant recoils effects in the environment. The characterization of the domain of 
validity of the stochastic treatment is a complicated question which requires a detailed 
knowledge of the environment dynamics [P^]. In the case of metric fluctuations, this is of 
course beyond our reach. But the crucial point is that, if there is an intermediate regime in 
which the particle weakly interacts with a large number of modes, it is sufficient to work 
with an appropriate simplifled model. In our case, universality comes essentially from 
the exponentially growing Doppler effect encountered in backward propagation: how ever 
small is the critical length it will be reached in a logarithmicly short (advanced) time. 
The only condition it must obey is to be much smaller than the Schwarzschild radius. 

The paper is organized as follow. In Section 2, we describe our fluctuating metric 
and we study how outgoing null rays propagate in it. The metric we choose results from 
the collapse of a massive spherically symmetric null shell which is followed by a small 
additional oscillating and infalling null flux. 

In Section 3, we analyse the propagation of wave packets of a massless scalar fleld in 
an ensemble of fluctuating metrics. That is, we flrst obtain the propagation of a given 
initial wave packet in each realization of the geometry and then define the mean wave 
packet by averaging over the fiuctuating part of the metric. Since the linearity of the 
problem is maintained, we formulate the problem as in a S-matrix language and show 
that this matrix is diagonal in Fourier components. From this equation one immediately 
sees that only the high frequency part of the spectrum is affected by the fluctuations. 

In Section 4, we flrst analyse the Green function in the initial vacuum and show 
that the metric fluctuations do not signiflcantly affect its short distance behavior thereby 



^ There are several works in the hterature wherein stochastic ensembles of black-hole metrics have been 
considered, see e.g. 0. However, to our knowledge, none of them confronts the role of ultra- high 
frequencies occuring in Hawking radiation. 
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guaranteeing unmodified properties of Hawking radiation. Tliis is verified by computing 
directly the asymptotic flux of energy of quantum radiation. 

In Section 5, we compare the backward propagations from of a wave function 
describing a typical Hawking quantum in two cases: in our fluctuating metric and by 



modifying the dispersion relation as done in [p!5|, |I6|. The similarities are manisfest. We 



conclude the paper by making comments on the apparent breakdown of Lorentz invariance 
which appears in these effective theories. 



2 Fluctuating Black-Hole Geometry 
2.1 Metric ansatz 

In this article we shall only consider spherical modes of metric fluctuations propagating 
in a spherically symmetric background. The most general spherical metric can be written 
in the form 

_ ^^A ^^B + ^^2 ^ ^2.1) 

where A,B = 0, 1, and 'Jab and R are functions of x^. The metric of a black hole of mass 
M formed by the collapse at f = of a massive null shell with mass M is 

dS^ = {AMfds^ , ds^ = -A dv^ + 2dv dr + r^d^l , (2.2) 

where in the absence of fluctuations 

A = Ao{r,v) = l~^, (2.3) 
2r 

{}{v) being the Heaviside step function equal to 1 for positive argument. For further 
convenience we have introduced the dimensionless coordinates {v,r), so that R = AMr 
and 4Mv are the radius and the advanced time in units where G = c = 1. A conformal 
Penrose- Carter diagram of the spacetime is shown in Figure 1. 

The most general metric perturbation preserving the form ( p.l|) of the metric is de- 
scribed by four functions of x^: 6r and S'Jab- The remaining coordinate gauge freedom is 
generated by infinitesimal coordinate transformations ^^{x). We fix the gauge by putting 

5r = 0, 5-frr = . (2.4) 

The perturbed metric can be written in the form 

ds"^ = '^{-Adv^ + 2dvdr) + r^dcul, (2.5) 

with 

^ = 1 + 5^, A = Ao + 6A. (2.6) 

Upon restricting attention to the propagation of radial null rays, the 2-dimensional 
conformal factor \1/ plays no role. In the following sections, we shall study the propagation 
of s- waves in the fluctuating black hole geometry. We shall demonstrate that, for s-modes, 
\1/ also drops out of the 4-dimensional Dalembertian. Hence only the function A will be 
relevant for us. 
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Figure 1: Conformal Penrose-Carter diagram of tlie unperturbed collapsing black hole geometry. 
The massive (of mass M) null shell is represented by a heavy solid line at u = 0. 



To further simplify the problem, we assume that the metric fluctuation SA is composed 
only of infalling radial null modes. Thus it is of the form 

6A=-^mKv)- (2.7) 
so that the perturbed metric is given by ( |2.5| ) with 

^ = ^„ + ,^=l-!W±iM. (2,8) 

2r 

For \1/ = 1 this is a Vaidya metric. The function encodes the light-like infalling 
fluctuations. As in eq. (f^.3D, the step function in relation ( ^.81 ) indicates that the black 
hole results from the gravitational collapse at f = of a massive null shell with mass 
M, and that there are no fluctuations prior to the collapse of the null shell. Therefore 
spacetime is flat to the past of the null shell. 



2.2 Stochastic variables 

To introduce the stochastic variables in simple terms, we postulate that fi{v) possesses a 
discrete^ and non-degenerate Fourier decomposition: 

) = ^ sm{ujv) + cos{tuv)] = ^ /Iq sm{ujv + 0^^) . (2.9) 

The equality is obtained by using polar coordinates in the (/ii, /X2)— plane 

/^o = \/(/if)' + (/i2 )' , tan0^ = 4- (2.10) 
/^i 

^ A discrete spectrum arises for example in York's approach |^ based on the quasi-normal modes of 

the black-hole metric. However, the spectrum due to other fields can be continuous. The results of the 

present paper can be easily adopted to this case. It is sufficient to replace the discrete sum, by an 

integral, J (iu;z/(u;), where v{uj) is the number density of fluctuation modes. 
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We assume that the (real) amphtudes and ^2 stochastic variables taking range 
from —00 to 00. For simplicity we postulate that they are stochastically independent 
and that, for any frequency oj, the dispersions of ii"^ and /ig are equal. Thus, there is no 
preferred value of the phase (p^ in the (/x'J', n^) plane. In this case, p^^fi^), the distribution 
function for the amplitude /Xq, satisfies the normalization condition 

2tt poo 

dcj)^ / dfx'^ fx'^ pUfJ^"^) = 1. (2.11) 
Jo 

Later in the text, we shall assume that the amplitude Hq is a Gaussian variable whose 
distribution function is equal to 



P-(^o) = TTZ^ exp 



2^2, 



f2.12) 



27ra2 

u. 

The coefficient determines the dispersion of the amplitude /Xg- 

<(m-)2>^= 2dl, (2.13) 

where O^^ represents the average calculated with the distribution Puj{Pq)- 

In what follows we shall have to deal with observables depending on the fluctuating 
geometry which obey the following factorization condition 

Q = \{QMAu). (2.14) 

For these observables, we can consider each sector labeled by uj separately. It is then 
useful! to introduce the successive averages: 

QM) = ^ r d(f)u. <5c.(/io, M , (2.15) 



27r Jo 

<Q^>^= 2tx / rf/i^ /i^ p^(/i^) Q^(/i^) , (2.16) 
Jo 

<Q>=n<Q->- ■ (2-17) 

The first equality gives Qa;, the average of over the stochastic phase (j)^^. The second 
one gives the result of averaging over the amplitude p^. Finally, the overall ensemble 
average of the observable Q, is given by (|2.17|) . The order in this averaging procedure 
follows from the fact that to perform the first average, one simply has to assume that there 
is no prefered direction in 0^. For the second instead, we need to know the distribution 
Pj^. And for the third one, we must know the whole spectrum. 

To determine the actual physical spectrum of metric fiuctuations around a black hole 
horizon is an extremely complicated problem requiring a theory of quantum gravity, see 
P, |18| for attemps to characterize this spectrum. In the present paper we shall not use 
any specific form of the spectrum. As we already mentioned, the effects of the metric 
fiuctuations hardly depend on its exact form. That this is the case will appear progres- 
sively in the paper. Our only assumption is that the amplitute of metric fluctuations are 
much smaller than the gravitational radius of the black hole. In our dimensionless units, 
this gives a <^ 1. This should be true for black holes of mass M much greater than the 
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Planckin mass nipi (in 0, the estimate dimensionless amplitude scales as a ~ (rripi/M), 
whereas in a ~ (mpi/M)^/^). 

Because of the factorizability of the operators and the statistical independence of the 
amplitudes, it will be sufficient to consider only a single fluctuation mode. To simplify 
the notations we shall drop the index u in the amplitude and in the phase. That is we 
shall work in the fluctuating geometry (p.2|), (|2.8|) with 



/i(f ) = fiQ sin(ti;f + (p) , (2-18) 



with /io <^ 1. We call the metric ( |2.2| ), ( |2.8| ) with fi{v) given by ( |2.18|) a realization of the 



fluctuating geometry. By averaging over (j) and /xq we thus assume that we are dealing 
with an ensemble of such realizations. 

Finally, it should be stressed that 0, the phase of the fluctuations, has physical mean- 
ing. Indeed, the gravitational collapse of the null shell singles out a special moment of 
time (e.g. the moment of formation of the horizon), and fluctuations with different phases 
are not equivalent. Therefore is a stochastic variable and when computing the statistical 
average, integration over it is to be done. 



2.3 Null ray propagation in a fluctuating geometry 

To characterize the propagation of radial null rays in a given realization of the geometry 
we focus on the relation 

V = V^{u) (2.19) 

between the moment v of advanced time when the null ray was emitted from J'~ and the 
retarded time u when it reaches J'^. The above equation contains a simplifled notation 
since V^(m) also depends on the amplitude /xq and the frequency u - this is evident in eq. 
(g) below. 

In the late time regime, i.e. m ^ 1, we have (for details, see 



vAu) 



1 + 



sm 



+ 00 + 



1 + ^ cos( 



+ 



UJ 



)) + 



: sm 



+ (po + uju - V5r(t^)) 



+ 0(/xo'), (2.20) 



where 0o = arctana; and where (fri^^) and q{uj) only depend on the frequency of fluc- 
tuations. In this Section and the next one, we shall work with a simplifled version of 
V(i){u): Since the last two terms inside the large parenthesis in ( |2.20| ) are much smaller 
than 1 [fiQ <^ 1), we shall simply omit them. In Section 4 instead, we shall take them 
into account and also include all quadratic terms in fiQ. 

The simplifled version of relations ( p.l9| ) and ( |2.2CI|) is of the form 



w 



WJu) 



(2.21) 



W^{u) = Wo sin(0 + 0o) + e" 
where we have introduced for later convenience 



w 



-1-v. 



Wo 



[2.22) 



(2.23) 
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For a given realization of the geometry and to first order in /xq fhe event horizon is 
given by the equation (see eq. (3.10) in 0) 

rTi^) = ^[1 + ^0 sin(a;t; + </> + </>o)] • (2.24) 
It crosses the collapsing null shell, f = Oorw = — 1, at the radius 

'^r(O) = ^[l + ^osin(0 + 0o)]. (2.25) 

Being traced backward in time the null geodesic giving rise to the horizon enters the flat 
spacetime region inside the collapsing shell, bounces at r = 0, and finally reaches 
with w (defined by (|2.23|) ) lying in the domain 

wei-wo,wo). (2.26) 

Therefore, for all (p, radial null rays emitted from J'~ at advanced time w < —wq fall into 
the singularity r = 0. On the contrary, radial null rays emitted with w > wq always reach 
The time of their arrival to J'^ lies in the interval 

mG(u_,u+), (2.27) 

where u± are 

u± = -\n{wTWo) . (2.28) 

Finally, the rays emitted in the interval ~Wq < w < Wq reach J'^ only for some values of 
(f). Since for w in this interval there always exists a phase such that the ray propagates 
along the horizon, the moment of arrival at J'~^ varies from u_ = — ln{w + Wo) to u+ = oo. 

Figure 2 shows the behavior of radial null rays in the fluctuating black hole geometry 
characterized by a single mode of frequency u. Ray a represents the event horizon, while 
rays b and c are rays which reach J'~^ or remain trapped inside the horizon, respectively. 
When they are traced back in time all these rays pass through the center of symmetry 
r = 0, and go to JT"" along lines v = const. It should be noted that the conformal 
Penrose-Carter diagrams for any given realization of the fluctuating black hole geometry 
is similar to the one shown in Figure 1 even though the event horizon no longer obeys the 



equation r = 1/2, see eq. (|2.24|) . 



3 Wave Propagation in a Fluctuating Geometry 
3.1 ^-pulse propagation 

3.1.1 Wave propagation in a given realization of the geometry 

In this article, we study only s-modes of a minimally coupled massless scalar field x- We 
introduce as usual (f = rx- Then, the four-dimensional Dalembertian equation = 



when computed in our fluctuating metric (|2^ ) gives, see e.g p3 



(2)n _ ^ 



ifi = (3.1) 
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Figure 2: Null rays propagation in a fluctuating geometry of the black hole. 

Hence \l/ defined in eq. ( |2.5| ) plays no role. Moreover, upon neglecting the centrifugal 
quantum potential, one obtains the equation 

(2)n^ = 2d^dr(fi + driAdrip) = . (3.2) 

When adopting this equation, one works in the geometrical optics approximation. There- 
fore, the spherically symmetric perturbations of the metric affect the global properties of 
the solutions of ( p. 21 ) only through the glueing of the null characteristics encoded in Vif>{u) 
given in eq. (|2.22| ). By "global properties" we mean properties which can be measured 
on J'~^. We shall not compute the local value of the field near the event horizon in the 
fluctuating geometry. This would require the knowledge of the local description of out- 
going null geodesies u^{v,r) for arbitrary v and r. Even though we restrict ourselves to 
global aspects, we shall see that the notion of a smeared horizon naturaly emerges. This 
should cause no surprise since the very definition of event horizon is global. 
We shall denote the value of the solutions of ( p.2| ) on by a capital letter 

$^ = y^\j± . 

Using the coordinate w defined in ( |2.23| ) we call $^(w) the initial value (or image) of the 
solution (f on J'~ , and $"*"(«) the final value (or image) of (f on J'~^. For a fixed geometry, 
i.e. for a fixed /iq and 0, the knowledge of ^~{w) uniquely determines its image on J'^. 
However, for backward propagation from J'~^ to J'~ , this is not true since J'~^ is not a 
complete Cauchy surface. 

Let us consider what happens to a wave packet propagating in the black-hole fluctu- 
ating geometry. The wave packet is sent from J'~ where it has image $~(w). Because of 
the linearity of the problem, it is sufficient to know the evolution of the following 5-like 
pulse sent from J'~ 

A-{w\w') = 5{w -w') . (3.3) 



9 



We call A^(m|w') the image of this pulse on in a given realization characterized by 
the value of 0. Since we work in the geometrical optics approximation, it is equal to 

AX{u\w')=b{WM-w'), (3.4) 

where W^[u) is the value of w at which a radial null ray has to be sent from J~ in order 
to reach at retarded time u. 

In terms of A J, the image of the wave packet on in this realization is given by 

/CO 
dw^-[w)A-l^u\w). (3.5) 
-co 



3.1.2 Averaging over the phase 

Using the notation introduced in eq. ( ^.151 ) we denote by A+(m|w) the average of A^(m|w) 
over the phase. Throughout the paper the "bar" over a quantity indicates that this average 
has been taken. In terms of A"*", the average image of the wave packet on ^ is given by 



$+(«)=/ dw^-{w)A^{u\w). (3.6) 

J —oo 

Thus our primary goal is to calculate I\'^[u\w). By definition it is given by 

A+(m|w) = -^ r b{W^{u)-w)d<^. (3.7) 

ZTT JO 

We remind the reader that this equation follows from the hypothesis that all values of 
are equally probable, see section 2.2. To calculate this integral, we use the following 
property of 5-function: 

rrf0<5(F(0)) = ^-i— , (3.8) 

where 0^ are the roots of the equation -F(0) = lying in the interval (0, 27r). That is, we 
simply need to find those solutions of the equation W^{u) = w which lie in the interval 
(0,27r). 

Using the simplified expression ( |2.22| ) we obtain^ 

A+ M w = \° ^ ' . 3.9 



t^\IWq — {e — w) 

For a fixed w, expression (|3.9| ) determines A+(u|iy) for u G (m_,m+), where 

M_ = — ln(w + Wo) , (3.10) 
i-Hw~wo), for w>wo; 

oo , lOT — Wo < W < Wo . 



■^Notice that A+ is highly non Unear in wq. However it results from eq. ( 2.22 ) which has been 
linearized and simplified. Therefore, one must question the validity of its non- linear depend ence in wq. 
This question is addressed in Appendix A. The outcome of the analysis is that one can trust ( [3.S| ) in the 
non-linear regime since it correctly gives the dominant non-linear effects. The fact that the higher order 
terms in /hq of are all irrelevant indicates that there is an underlying universality of our results. 
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For u ^ A^{u\w) vanishes since no G (0, 27r) satisfies W^{u) = w. Similarly, 

for a fixed value of u, A'^{u\w) does not vanish only for w G {—wq + e~", wo + e~"). 

In a similar way we can consider propagation backward in time. It is again sufficient 
to make the calculations for a 5-like pulse. We call 

+ A{u\u') = 6{u - u') (3.12) 

the image of the pulse on J''^. Then, for a given realization of the geometry it has the 
following image on J'~ 

-A^{w\u')=5iu'-U^iw)). (3.13) 

Using eq. ( p.22| ), one has 

U^{w) = -\n[w - wosm{(f) + . (3.14) 
By averaging ^A^ over we get 

-A{w\u) = e-"A+(M|w) , (3.15) 



where A^{u\w) is given in ( |3.9| ). The only difference is the Jacobian e~" which relates the 
delta functions defined on J'~ in terms of w and on J'~^ in terms of u. Then, the image 
on J'~ which results from the backward propagation of a wavepacket which has the image 
on J'~^ is given by 

/oo _ 
du<t>+{u)-A{w\u) . (3.16) 
-oo 

It should be noted that in order to get this relation we assumed that, for any realization 
of the geometry, no signal was propagating backward from black hole interior. This point 
will be further discussed in section 3.3. 



3.1.3 Averaging over the amplitude 

The average value of an observable over an ensemble of amplitudes of metric fiuctuation 
IIq characterized by the distribution Pi^(/io) is given by eq. ( p.l6| ). In order to compute 
the amplitude average of A"*", we need to know p^j. In the case of a Gaussian distribution 
( p.l2|) , the average can be easily performed. Since the impact of the metric fiuctuation 
on the field is expressed by wq = po/V^ + uj"^ instead of the origin amplitude /io, we 
introduce the new distribution function 

1 w"^ 

Pc^iwo) = - — -exp{-—^) , (3.17) 



where 



27ra2 " ^ 2(t2 



a^ = -^^. (3.18) 



It is easy to check that 

i-OO 

27r / dwoWop^{wo) = 1 , < Wq >^= 2ol . (3.19) 
where the average is now taken over Wq with the distribution function Puj{wq). 
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Using using ( p.9| ) and ( p.lTj ), the average of A"*" gives 



< A+(m|w)>, 



cr„,V27r 



exp 



2^2 



(3.20) 



Thus the Gaussian character of p^j, the distribution of metric fluctuations, is preserved 
and now characterizes the distribution of hght rays. The same result applies to backward 
propagation since one has <~A{w\u)>^= e"" <A~^{u,w)>^. 

At fixed u, < gives the image on J'^ of the given -u-ray sent from J'^. This 

image gives the probability that this ray reaches the interval dw centered on w. For 
u = oo, < A+ therefore gives the image of the event horizon. Since we get a 
Gaussian distribution centered at w = 0, the position corresponding to the horizon in 
the unperturbed geometry, this indicates that the horizon is smeared by the stochastic 
fluctuations. And eq. (|3.2CI|) shows that the global properties of the field propagation are 
sensitive to this smearing of the event horizon. 



3.2 Scattering operator in a fluctuating black-hole geometry 

The relations (|3.4|), ( p.9| ) and ( |3.20| ) solve the classical scattering problem in our fluctu- 
ating geometry. We now demonstrate how this problem can be solved in a more formal 
way. This more formal approach will be very useful for understanding the modifications 
of the propagation induced by the metric fluctuations. 

In order to simplify the algebra it is most convenient to introduce a new coordinate 
on J+ 

y = e-^. (3.21) 
The reason is the following. In terms of y, eq. ( p.22| ) becomes 



w = W^{y) = y + WosH<p + (t)o). (3.22) 

Therefore, the effect of the metric fluctuations is simply to shift y with respect to w. The 
image of a wavepacket on can be considered as a function of u or of y. In order 
to avoid confusion, we shall keep the notation $"''(«) for the function of u and shall use 
the notation ^^{y) whenever the image is considered as a function of y. 
In terms of A^{u\w) given by (|3.9| ) takes the form 



A-(,H = i^^p^iE^, (3.23) 



and relation (13.61) becomes 



- , 1 f^o+y dw $ iw) 

^Hy) = Z I ^ I , / (3-24) 

-wo+y y'ty^ - [y - wY 



IT 



According to our definition (|3.21|) the coordinate y takes only positive values. But we 



can use relation ( |3.24 ) to define ^~^{y) (at least formally) as a function on the complete 



y axisQ. In this case, eq. ( |3.24| ) defines a linear operator acting on the space of functions 



defined on the real axis. We call this operator D. Thus ( ^.241 ) takes the form 



= (3.25) 



*The geometrical meaning of negative values of y will be presented in the next subsection. 
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A formal representation for this operator can be easily obtained by exploiting the fact 
that eq. ( ^2^ is a linear map. Indeed, for any realization of the geometry we have 



= + «7o sin(0 + 0o)) = e"'°^^^(<^+^»)^ . 



(3.26) 



Here d is the operator of differentiation with respect of the argument of the function. 
This relation can be rewritten as 



where 



^WQ sin(0+</>o)9 



(3.27) 
(3.28) 



Notice that this shift operator bears some similarities with that introduced in 

Let us first perform the average over the phase. Using the integral representation of 
the Bessel function of zero index 



Mb) = - 

TT Jo 



we get 



D 



2tt 



27r 



D</, = Jo{-iwod) 



(3.29) 



(3.30) 



Let us now show that ( p.23|) is a direct consequence of ( p. 301 ). For this purpose it 
is convenient to adopt the Dirac notations and to write the functions and $~ as 
ket- vectors |^^) and |$~), respectively. Then, in the "coordinate" representation, we 



have 



and equation ( |3.25| ) takes the form 

$+(x) = 



dx'{x\D\x) ^"(x') . 



(3.31) 



(3.32) 



Relation (|3.30| ) implies that the operator D is diagonal in the "p-representation" . This 
directly follows from the linearity of eq. ( |3.22| ). Explicitly one has 



{p\D\p') = 5{p-p') Joiwop) . 



Then, using 



we get 



{p\x) 



exp (ipx) 



(3.33) 
(3.34) 



dp dp'{x\p){p\D\p'){p'\x') = — dpe-'^''-^^P Jo{wqp). (3.35) 

-OO J —CO ZTl J —OO 

Calculating the integral (see e.g. [^) we get 

1 ^{wl -{x- x'f) 



(x|D|x') 



(3.36) 
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Thus {y\D\w) equals A+(?/|w) given by ( |3.23| ). 

Up to now we have only been dealing with the stochasticity connected with the phase 
0. Let us discuss what happens when we average the operator D over the amplitude of 
the metric fluctuations. We again assume that the probability distribution is Gaussian 
( p. 171) . Using the relation (see ||26| p. 393, eq. 13.3.1) 



°° " 1 .,2 



dzze-' Jo{az) = -e"" , (3.37) 
2 



we get that the average of D over the fluctuation amplitude is 



<D>^=exp ) . (3.38) 



This simple form results from the linearity of the map ( p.22|) and the Gaussian character 
of <A+{u\w)>^ (PI). 

The generalization of these results to a spectrum of metric fluctuations with different 
frequencies is straightforward. Indeed, since < D has a simple exponential form, it 
enjoys the factorization property ( |2.11| ). Thus, using ( |2.17| ), its total ensemble average 



<D>=exp , (3.39) 



where 



cr 



cff 



^ - E ^^ (3-40) 



Eq. (|3.40D shows that the effect of the whole spectrum of metric fluctuations is to give 



rise to a single length which weights higher order derivative terms. This shows that the 
details of the fluctuations spectrum play no significant role for the scattering operator 
<D>. 

In resume, finding the image on J'^ of a wavepacket determined by its data $~ 
on J'~ is a "classical" S-matrix problem. In the absence of metric fluctuations, the linear 
operator So relating $~ and is trivial (equal to 1) if we use the coordinates y and w. 
In the presence of a given realization of the metric the operator is given by eq. ( p.28| ). 
When one computes the "mean" propagation by considering a stochastic ensemble of 
metric fluctuations, this simple relation is modified and takes the form ( p. 301 ), ( p. 381 ) or 



( p.39|) according to the ensemble of metric fluctuations that one considers. 

The extremely simple form of the operators D in the "p" -representation allows one to 
make a few general observations. In particular, we have 

<$+(p)>= {p\ <$+>) = e-"offPV2$-(p), ^3 4^^) 

Thus only the high frequency components (i.e. of the order of a~g^ and greater) of the 
initial wave packet are strongly affected by the fluctuations of the geometry. Therefore 
the forward propagation of any smooth (i.e. of Fourier content much below cr~J^) wave 
packet defined on J'~ will not be significantly affected by the metric fluctuations. In 
other words, for classical black hole physics, the metric fluctuations are irrelevant if, as 



indicated in |I8[, their "mean" amplitude a^g is of the order of the Planck length or 
smaller than it. 
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3.3 Backward in time scattering 



Eventhough eq. (|3.22|) is perfectly symmetric in w and y, backward propagation is dra- 
matically affected by the metric fluctuations. 

To settle the discussion, we first clarify the geometrical meaning of negative y. To ease 
this analysis, we start with backward propagation in the absence of fluctuations. In this 



case, for positive y, eq. (p.22|) gives w = y. However, since J'^ is not a complete Cauchy 



surface, we need to consider the union of J^~^ and the whole event horizon m = oo in 
order to have a complete Cauchy surface. Thus we must introduce a coordinate along the 
horizon. The simplest choice consists in considering the negative values of y defined again 
hj y = w. Indeed the negative y axis so defined covers the horizon from r = 0,w = 
inside the collapsing shell till w = — oo. Thus, the real axis y G (—00, 00) forms a complete 
Cauchy surface and the functions determine their image on J'~ for all w. 

This procedure also applies for any given realization of the fluctuating geometry. In- 
deed, the event horizon, when continued backward for negative v in the inside flat geom- 
etry ( |2.3| ), reaches r = at advanced w time ^^(y = 0) = Wq sin(0 + 0o)- Therefore the 
negative half line y G (—00, 0) defined by eq. (|3.22| ) still covers the whole horizon and 
y G {—00, 00) forms a complete final Cauchy surface. Since this is valid for any realization, 
it is meaningful to use the coordinate y G {—00, 00) after having averaged over 0. 

For y > 0, ^~^{y) gives the value on J'^, while for y < it gives the value on the 
horizon. For regular ^~^{y), in virtue of the symmetrical role played by y and w in eq. 



( p. 221 ), the averaged image on J' is determined by the same scattering operator D which 
governed forward propagation. In the case of the full ensemble average, one has 



In Fourier transform with respect to w and y this gives 

<<l>"(p)>= exp 



(3.42) 



(3.43) 



Of special interest for studying Hawking radiation are the final images such that no 
incoming field emerges from the horizon for any realization of the geometry. They are of 
the simple form 

^tutiy) = m^''iy)- (3.44) 

Unless ^~^{y) vanish sufficiently rapidly when y 0, these modes are singular aX y = 0. 

This problem will be studied in detail in Section 5. It reveals the important role 
played by the fluctuating horizon geometry for backward propagation. The asymmetry 
between backward and forward propagation comes from the fact that the inertial time on 
J'~^ which characterizes out-frequencies X = du is u and not y = e~^. Then, the so defined 
out-frequencies are exponentially blue-shifted when propagated backward near the event 
horizon. This purely kinematical effect is at the origin of the trans-Planckian problem 
and has here dramatic consequences since higher derivative terms are present. Indeed, 
eq. (|3.43|) when applied to out-functions (|3.44|) which vanish for negative y gives, in the 
position representation. 



(w)>= 



exp 



- In w 



(3.45) 



The dramatic consequences can now be seen: how ever smooth is the final data $"'"(«), the 
fluctuations of the geometry will inevitably affect its backward propagation if it is centered 
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around a sufficiently late retarded time. Moreover if it does not vanish sufficiently fast 
(i.e. faster tliari0 e~") when one approaches the horizon, the asymptotic behaviour of 
^ D ^ intervenes. 

We shall further analyze these points in quantum mechanical terms in Section 5. 
Before presenting this material we shall show that the asymptotic properties of Hawking 
radiation are not significantly affected by the fluctuations when their average amplitude 
CTefi is much smaller than the Schwarzschild radius (= 1/2 in our units). 

4 Hawking Radiation 

4.1 Green function in the m- vacuum 

The simplest way to understand why the metric fluctuations do not significantly modify 
the asymptotic properties of Hawking radiation consists in analysing the Green function 
(more specifically, the positive frequency Wightman function ||2^) evaluated in the initial 
vacuum state. Indeed, as shown in when the short distance expansion of this function 
evaluated near the event horizon reduces to the standard (Hadamard) behavior. Hawking 
radiation obtains on . 

In the absence of fluctuation, the m-vacuum is the vacuum state with respect to 
positive frequencies defined on J~ in terms of the inertial advanced time derivative id^. 
Therefore, for late time, the w-part of the Green function evaluated in this vacuum is 
controled by V[u) = — 1 — e~". Explicitely, one has 



The thermal and steady character of Hawking radiation follows from the exponential 
relation between V and m, which remains unchanged for arbitrary large time m, and from 
the fact that for small V intervals [5V <^ 1), the in-Green function behaves like ln(5V^. 
Therefore these are the two conditions which should be met in order to see if the metric 
fluctuations lead to modifications of Hawking radiation. 

In a given realisation of the metric fluctuation, one simply replaces the unperturbed 
relation V{u) = — 1 — e~" by its modified version eq. ( p.20| ). Upon considering the en- 
semble of geometries, one averages the Green function over the ensemble. After averaging 
over phases, the mean function is thus given by 



^This condition on the decrease of the wave packet already appeared in the literature P9| in the 
Unruh detector context (see also for its transposition in black hole physics). In the Unruh case, if 
one imposes that the fluxes of energy emitted by the accelerated detector be regular in an inertial frame, 
the coupling between this detector and the field must decrease faster than e"""^ where a is the acceleration 
and r the detector's proper time. 




]- \^{V{u) - V{v!) + ie) = ln(-e-" + e""' + %e) . (4.1) 
Ivr 47r 




(4.2) 
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First, we notice that when using the simphfied expression ( |2.22| ) Hawking radiation is not 
at all affected since the constant shift cancels out. The reason is clear: in each realization, 
two nearby points are shifted by almost the same amount, see Figure 2. 

Therefore only higher order effects might affect Hawking radiation. To compute them, 
it is again usefuU to use y = exp(— m). In terms of y, using eq. ( |2.20| ), one has 



where 



V^{u) - V^{u') = F^{y') - F^{y) + 0{wl) 



F^{y) = y{l + acos(0 + 20o) + 6sin(0 + - ^^Iny)) • 



(4.3) 



(4.4) 



Here ipo = (po " y^r- The explicit form of the coefficients a and b can be easily obtained 
from eq. (|2.20|) . They are not important for us now. It is sufficient to mention that both 
these coefficients are of first order in wq. For close points y' = x + z/2 and y = x — z/2 
(z<^l) one has 



F^{y') - F^{y) 



1 + H Cn(0)^ 



2n 



n=0 



(4.5) 



All the coefficients c„ are of the first order in wq since F^ has been linearized. In particular, 
the first coefficient 



co(0) = ci cos(0 + 20o) + b sin(0 + ip^ — ojlnx) — hu cos(0 + ipo — uj\n.x) . (4.6) 

is a linear superposition of sine and cosine whose arguments are linear in 0. It is easy to 
check that this is also true for all other coefficients c„(0). Then, using ( [4.3|) and ([4.5|) , we 
can expand ln(V^('u) — V^{u') + ie) in powers of wq and conclude that after averaging over 
one has 

^^"(m, u') = G^"(m, u') + wlHiu, u') , (4.7) 

since the averaged value of all first order terms in wq vanishes. The crucial point is that 
H{u,u') is finite when u — >■ u'. This implies that the corrections to Hawking radiation 
are at least of second order in wo and non-diverging in the late time regime. Upon 
considering the average over the fluctuation amplitudes the same result holds. These 
points are explicitely verified in the next section. 



4.2 Quadratic corrections to Hawking Flux 



The (quantum mechanical) mean energy flux of Hawking radiation measured at J7^, 
dE/du, was obtained in ^ in the case of the oscillating Vaidya ( |2.8| ) metric with the 
phase (f) in eq. (|2.18|) put equal to zero. The corrections were computed up to second 
order in the fluctuation amplitude /iq by using the 2D model based on eq. (|3.2|). In that 
case, the energy flux can be decomposed into the sum of a permanent part {dE/duY""'^ 
and a fluctuating part {dE / duY'^"^ 



dE/du = (dE/duY^''^ + (dE/du) 



The permanent part is equal to 

{dE/duY"''^ = 



K 



487r 



I + -wl{l + uj^)q\uj) + wl 



(4i 



(4.9) 
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where k = (SvrM) ^ is the surface gravity of the unperturbed black hole, and 



q{uj) 



- 1) 



(4.10) 



In the absence of metric fluctuations, the permanent part of the flux reduces to its usual 
value, /t^/487r. The corrections due to these fluctuations are second order in w^. The 
fluctuating part {dE / duY"^"^ is a linear combination of terms oscillating with frequency uj 
and 2uj and, by definition, vanishes when integrated over the retarded time u. 

We shall now compute this flux in a stochastic metric by averaging its value over the 
phase and amplitude of the metric fluctuations. We first study the influence of a single 
stochastic phase 0. In particular we shall check that averaging over this phase and over 
u does not affect the value of permanent part of the flux. For this purpose, as in 0, we 
calculate dE/du from the 2D expression 



dE_. 
du 



127r 



{dV^/duy/'^ [{dV^/du)-'/' 



(4.ii; 



where Vtf){u) now depends on the phase 0, see ( p.20| ). Since we want to compute the 
quadratic corrections in wq, we need to know V^{u) up to second order in wq. In the late 
time regime and to this order in wq, V^{u) can be decomposed according to the dependence 
of the various terms in u. Using eqs. (4.51-60) in 0, one gets 



VJu) 



Vo(0) 



Wn 



Vr{<j)) + V2{u,<j))-^u 



(4.12) 



By definition Vq and Vi do not depend on u. Note that Vq, Vi and V2 parametrically 
depend on the frequency u of the fluctuations. 

For a given 0, i.e. for a given realization of the geometry, one easily sees from (|4.11| ) 
that dE/du does not depend on the values of Vq and Vi. Therefore one is free to put 
Vq = and Vi = 1 and to work with V^{u) as 



VJu) 



Wn 



1 + V2(U,(P)--^U 



(4.13) 



The novelty with respect to eq. ( p.20|) arises from the fact that V^{u) is now developed 
up to second order in wq. As V2{u, (p) has no 0-th order term in wq it can be decomposed 
as 

V2{u, 0) = WqV2,i{u, 0) + wlV2,2iu, 0) . (4.14) 
The first term is already known, see eq. ( |2.20| ), 



V2,i(m, 0) = g(t^) sin(0 + 00 + - v^r) , 



(4.15) 



where V5r(i^) is a phase independent of u. V2,2(^, 0) is a complicated expression containing 
oscillating terms with respect to u and 0. It is of little interest to write it explicitely. In 
what follows we shall only need its averaged value with respect to since we are interested 
by the expression of the averaged energy flux to second order in wq. 
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In terms of these two functions, the energy flux reads 



dE 
du 



u,(p) 



K 



487r 



(4.16) 



where ' stands for d/du. 

Using as in Section 2 a "bar" to denote averaging over one gets 



dE 
du 



[u] 



487r 



(4.17) 



The first Wq term was already present in eq. ( [4.9|) and is due to the last term in ([4.12|) . 
The second one comes from quadratic terms in V2,i given by ( |4.15| ). The last term depends 
on u and is equal to 



Vi2 



2,2 



(1 — 2uo^)q{ijj) cos{uJU — (pr) + uj{3 — Auj) sin(c<j'U — (fr) ■ 



(4.18) 



These results show that after taking the averages over and over u one gets back the 
same permanent contribution to the Hawking mean flux as in a given realisation of the 
geometry, see eq. ( [4.81 ). The reason is that averaging over u erases all dependence in (p. 

Eq. ( [4.17] ) can be easily extended to a spectrum of metric fluctuations according to 
the lines presented in section 2.2. Limiting ourselves to the permanent part of ([4.17|) we 
obtain 

1 



dE' 
du , 



487r 



I + WqM ! + -(! + cj")g"(cu; 



(4.19) 



which gives the average over the phase of a single fluctuating mode of frequency uj. When 
using the Gaussian distribution p^, see ( |3.17|) , the average over the amplitude wq^u) of 
the fluctuating mode gives 



< 



'dEy 
du I 



4^[l + 2aS(l + l(l + .VM) 



(4.20) 



where ( p.l9| ) has been used. 

Upon considering the whole spectrum of fluctuations, we get 



'dEV^ 
du J 



487r 



1 + 25:^5 (l + 1(1+.; vh) 



(4.21) 



To perform this last summation requires the knowledge of the spectrum, here represented 
by the set of a^. 



5 Fluctuating Geometry and Trans- Planckian Prob- 
lem 

5.1 Hawking radiation for modified dispersion relations 

In this Section we shall establish the close analogies between the backward propagation 
in a fluctuating metric and the altered propagations which have been recently studied 
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and which result from the modifications of the dispersion relation in the high frequency 
regime. 

Before presenting the technical details, it is appropriate to recall the following points. 
These models have been introduced in order to show that the mutilation of the dispersion 
relation for frequencies higher than ujc-, which is the equivalent of in our case, in 
no way affect the (low energy) properties of Hawking radiation, namely stationarity and 



thermality. Following the original work of Unruh[p, |T5[, many models have been analysed. 
Their common property is that the Dalembertian is modified by the addition of higher 
derivative terms weighted by negative powers of uoc- In these works, the modifications 



have been inspired by hydrodynamics , electrodynamics in a dielectric medium p2 



field theory on a lattice theory []33|, string theory or by guessing what the physics near 



a horizon might be|jT6|. In all these models, the following properties obtain 

1. Fowardly propagated wave packets are unaffected by the modification of the dis- 
persion relation as long as their in-frequency content is much below the critical 
frequency ojc- 

2. No significant modifications of the asymptotic properties of Hawking radiation as 
long as the surface gravity satisfies k -C oJc- 

3. Dramatic modifications of backward propagated late-time wave packets of out- 
frequency A when the blue shifted value Ae" reaches Uc- 

It should be already clear to the reader that our effective propagation in a fluctuating 
metric possesses many similarities with these models. Our aim is now to establish the 
parallelism in simple and analytical terms. To this end we shall exploit the following facts. 

First we exploit the stationarity of the unperturbed background geometry by consid- 
ering the backward propagation of a plane wave of out-frequency A defined on J'^ . We 
shall consider its image on i; = rather than on J~ in order to ignore the backward 
propagation from f = till which is very much dependent on a model of a collapsing 
body and presumably irrelevant for black hole physics. The relationship with what we 
did in the former Section is straightforward since w, defined on corresponds to 2r — 1 
on f = in the dimensionless coordinates defined in Section 1 where k = 1. 

Secondly, we work in Fourier transform with respect to w as it provides an elegant 
characterization of in- vacuum. Since r is an afine parameter along v = const, —idr 
corresponds to positive frequencies measured by inertial observers when they cross the 
event horizon. Therefore, when infalling observers see no particles, it means that the state 
of the radiation field corresponds to vacuum with respect to frequency modes of positive 
—idr- Moreover, since dr\v=o = '^d^lj-, the initial vacuum with respect to w coincides 



with vacuum as seen by infalling observers, see [T^l • 

Lets now review how these concepts translate in mathematical terms and how they 
can be used to compute the consequences of modifying the dispersion relation at high 
frequencies. In this we shall present the "alternative" model of for its simplicity and 
its generahty. 



The 2D Dalembertian ( |3.2| ) for out-going modes with out-energy A propagating near 
the event horizon (r — 1/2 ^ 1/2) gives, in the unperturbed Vaydia metric 



;i - 2r)idr(px = 2id^Vx = 2\ipx ■ (5.22) 
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Along = 0, it is useful to express r in terms of w = 2r — 1. This leads to the simplified 
equation 

wdy^ipx = iX^x . (5.23) 

The general solution is of the form (px = A-t}{w)w^^ + B'd{—w){—wy^. The normalized 
out-mode describing the one particle state of energy A is given by 



fx 



w 



iX 



(5.24) 



On the other hand, the mode leaving on the other side of the horizon describes the partner 
of this Hawking quantum. 



In Fourier transform p = id^ eq. ( |5.23| ) becomes 



dp (pv^a) = -i>^Vx 



(5.25) 



In terms of p, the normalized in-mode contains only positive p and is thus given by 
_ ■(9(p)p~*'^~^/ \/ AttX. To obtain the Bogoliubov coefficients encoding the thermal fiux 
of outgoing quanta it suffices to inverse Fourier transform this mode and use the definition 
( ^.24|) . Explicitely one has 



dp 

y27r' 
Ti-^X) (6 



IW 



AX 



271 V4^ 



r(-zA) 



27T 



VAttX 



(5.26) 



where T{z) is the Euler gamma function and where e is small and positive. It fixes the 
relative weights of w^^ for positive and negative real valuesj^, as shown explicitely in the 
third line. This relative weight determines in turn the ratio of the Bogoliubov coefficients 
ax and Px- \ax/ Px\ = e^^- 

For later purpose, we notice that the same result can also be obtained the other way 
around, by Fourier transforming the out-mode (|5.24|) . In this case one gets 



°° dw 



e *P"'(^r*(w) 



V27r 
T{iX + l) {e + ip)-^^~^ 



y/2^ 

-i)V{iX + 1 



271 



i-p) 



-iA-l 



(5.27) 



The only subtility concerns again the prescription given by e. It arises this time from the 
fact that the final wave packet vanishes on the other side of the horizon. Upon Fourier 
transform, this fixes the relative weights on the positive and negative real p axis. This 
can be used again to determine the ratio of the Bogoliubov coefficients. 

We are now in position to modify the dispersion relation. In a fiat metric, this relation 
is simply p = X. To modify it, we write it as g{p) = X. The only condition that g must 
satisfy is that for small p is behaves as g{p) = p{l + 0{p/ujc)). Near the horizon, for 
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ujc ^ 1 and up to a normal ordering ambiguity of dp and g{p) which plays no role in a 
WKB approximation, eq. (|5.25|) becomes 

dp {9{p)^x) = -«A<^A • (5.28) 
The general solution of this modified equation is given by 

-iX f dp'/g{p') g-iA f dp'/g{p') 

(p^ = Ad{p) — + B'd{-p) — . (5.29) 

9VP) 9{P) 

In these terms, the three properties listed above are easily obtained. The first one 
follows from the definition of g{p) which deviates from linearity only for p > ujc- The sec- 
ond point is verified by using the fact that in-modes characterizing vacuum for infalling 
observers still contain only positive p. The proof goes as follows. The Bogoliubov coef- 
ficients are still determined by taking the inverse Fourier transform of (^\", given by eq. 



( [5.29[ ) wherein one puts A = l/y/inX, and B = 0. Then, one sees that for \w\ ^ l/ujc, the 



integral is dominated by values of p in the cis-Planckian domain p <C cOc- Irrespectively of 
the nature of corrections to the dispersion relation encoded in g{p), this locality implies 
that for these 'large' \w\, one recovers the un-modified propagation on each side of the 
horizon which characterizes the out modes. This in turn implies that one also recovers 
the usual Bogoluibov coefficients, for more details see |jl6[ . 



In order to discuss the third point one must choose the deviation from linearity in g{p). 
For our purpose, it is sufficient to consider g{p) of the form g{p) = pil+^p^ /ul). For ^ = 1 
one obtains sub- luminous propagation, for ,^ = —1 super-luminous propagation and for 
^ = i one gets dissipation without significant dispersion. These results are easily reached 
by constructing wave packets and analysing the locus where their phase is stationary. The 
main point is the following: when the phase of the function varies faster (slower) than the 
unmodified phase given by Alnw, one has sub-luminous (supra-luminous) propagation. 
In anticipation to what we shall get for metric fiuctuations, we say a few more words in 
the case of dissipation. In this case, the first order deviation for A ^ 1 is of the form 



^r(p) 



exp(-zA/^cip'^^ 



ip''^/u'i) 



V47rA 1 p(l + ip'^/ujl) 

xp^/2ul ^ 0{p^/uol)) (5.30) 



e 



This essentially corresponds to that will be obtained in the case of metric fiuctuations. 
The fact that dissipation (and not only dispersion) should physically occur for black hole 
was discussed in and the space time image of a packet built with waves of the type 
( ^.30 ) was schematically presented in Fig. 5, see also [^] for a super- luminous dispersion 



relation which leads to an effective dissipation of the wave packet. 

There is one more point which should be discussed. It concerns the physical relevance 
of considering propagation backward in time. A few remarks may illustrate its relevance^ 
First backward propagation provides the simplest tool to investigate the nature of the 
field configurations which give rise to Hawking radiation. In particular, it clearly estab- 
lishes that Hawking quanta emerge from trans-Planckian configurations when one uses, 

^However in the absence of a manageable theory of quantum gravity, we are still missing an explicite 
computation of backreaction effects which will settle this question in unambiguous terms. 
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as Hawking originaly did, the free Dalembertian for propagating the modes. Secondly, 
when considering S-matrix elements in a quantum dynamical framework, backward prop- 
agation always occurs since one also fixes the final state of the field, see [13, T^\. The 



simplest example is provided the Feynman in-out Green function. This function can be 
easily related to in modes characterized by an out-energy A. Indeed, one has 

^ ^ J <0,out\0,in> 

du e*'*'" Gin-outlw, f = C)o; w, = 0) (5.31) 

Because of its out-frequency content, the wave function is global in the sense that it 
entails the propagation from to the point where it is evaluated, here (w = 2r — 1, v = 
0). In particular when computed in a fluctuating geometry, this function encodes the 
effects of the metric fluctuations one encounters from to the point where it is evaluated. 



5.2 Backward propagation of owt-modes 

We now have all the tools to perform the comparison. To be as close as possible to 
what we just presented, we consider the image on f = of the monochromatic wave of 
out-frequency A. That is, the image on is = (Px'^iy = w) given in eq. ( p.24| ). 



Then we have to face the problem mentioned after eq. (|3.44|) : ^°x^{y) is not a smooth 
function but a distribution. This is of course a manisfestation of the trans-Planckian 
problem. For instance, its derivative with respect to y is ill defined on y = 0. Therefore, 
one must regularize ^p°^^^ before applying the scattering operator D on it. 

The simplest way to define the action of ^ D ^ on our owt-function is to work in the 
momentum conjugated to w. Indeed, the Fourier transform of ^p°^^^ is well defined in the 
high p regime, see eq. (|5.27|) . Then using eq. ( |3.43|) we simpy get 



« *,-.,(p) »= exp (-^ -^UV) = -P (-^) ^-^(^ + 'P)-"-' 



exp 



pVi\ ( ^)^+^) [e-V2^(p)p-^-i - e-V2^(_p)(_p)-a-i] . (5.32) 

Since the effect of the stochastic fluctuations is to multiply the wave function by an even 
function in p, the relative weight encoding the Bogoliubov coefficients is unaffected. This 
guarantees that the vacuum state with respect to p > leads to the usual properties of 
Hawking radiation, thereby proving point 2 above. Moreover, eq. ( p.32| ) confirms that 
the trans-Planckian problem is tamed: The high frequency content, i.e. the near horizon 
behaviour, is suppressed by a Gaussian factor, as for the dissipative case in the former 
Section, see eq. (|5.30|) . Finally, for p^a^g ^ 1, one recovers the usual expression of 
out- modes in terms of p given in ( |5.27| ), thereby proving point 1 above. 

To complete our analysis, we shall now determine the behaviour of <C ^ in 
spacetime. To this end, we inverse Fourier transform separately the two terms (positive 
and negative p) which appear in eq. ( p.32|) . Using eq. (1.11) in |Q, we obtain 



err 



« $(,)(«;) »= Z^{w/a^,) , (5.33) 
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(a) (b) 

Figure 3: Plot (a) gives the real part of Z\{x) as a function of x for A — 5.0. Plot (b) gives an 
higher resolution of the same function in the interval —2.0 < x < 2.0. 



'^^'^D,x{t - ix) - e-^^/^Dix{t + ix)] . (5.34) 



where 

2 smh(7rA) 

Here D^lz) is the parabohc cyhnder function. This function is related to the Whittaker 
function Wk,u{x) (see, e.g. p8|, p. 39) 

D,{z) = 2^/' (^^) W^./2+i/4,i/4 (y) . (5.35) 

In (|5.34| ) we have introduced e positive and infinitesimal in order to specify the phase of 
e — ix for positive and negative x. Figure |^ illustrates the behavior of Z\{x). We have 
plotted the real part of Z\{x) using Maple and relation ( ^.35| ). 



The function Di\{z) has the following asymptotic behavior for \z\ ^ |A|, see p7 



Aa(^)~^^V^'/'|i + ^^^ + 0(z-^)| , for-^<argz<| . (5.36) 

Using this relation, one verifies that for large negative values of x = w/a^g, Z\ vanishes. 
Instead, for large positive z behaves as 

^a(x) ~ o:^^ |i - + 0(x-^)| , (5.37) 

The limit x = w/a^f^ —>■ oo corresponds to the far from horizon region or to a^g 0. 
The last case corresponds to no metric fluctuations. In this regime (|5.371 ) reproduces the 
unperturbed out wave function given in eq. ( |5.24| ). The plot of the unperturbed out- wave 
function x*'^ is shown in Figure ^. By comparing Figures ^ and ^ one clearly sees that 
metric fluctuations strongly affect the back scattered wave for values of w close to smeared 
horizon. Indeed, close to the horizon, i.e. for -C A, the function Zx behaves as, see 

Zx{x) = (^i^^ |l + coth(7rA/2)e^'^/^xv^A - t/2 + Oix^X)^ (5.38) 
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Figure 4: Plot of the real part of the unperturbed out-wave x''*' for A = 5. 



The most important feature is the disappearance of the infinite trans-Planckian reservoir 
of fluctuations of equal amplitude which characterizes the unperturbed out-wave w"^^. 
This implies that backward propagation stops for any localized wave packet. In order 
words, the amplitude of the wave packet is dissipated in its backward motion once the 
wave packet enters the near horizon region a; < A. 

To show this, let us consider a rapidly modulated Gaussian wave packet which on J'~^ 
has the image 

= e-*^("-"°)e^^^i^ . (5.39) 

This image is localized around u = uq and has width b. In the limit b ^ oo one obtains 
a monochromatic out-plane wave of frequency A. Its Fourier decomposition is given by 

/) PCX) — 



^ ^ ' A/2vr ^ -00 

We first consider its backward propagation in the usual non-fiuctuating geometry. In this 
case, the plane wave e~*^" takes the form w^^ on f = 0. Thus the image on w = of the 
wave-packet ( |5.39| ) is 

$- (ty) = g*A(ln«>+?io)g-(liHi'+Mo)V(2&^) _ (5-41) 

For late time mq, it hugs the horizon = since it is centered around w = e~'^°. More- 
over it is highly blue-shifted since its mean frequency is equal to Ae"° for an infalling 
observer. Finally, no matter how large uq is, its maximum amplitude is still 1. There is 
no dissipation. This is the trans-Planckian problem. 

Let us now analyze how the fluctuating geometry changes this picture. In the presence 
of fluctuations, the image on f = of the plane wave e~*'^" is al^ Zx{w/ a ^ff) instead of 
^iA_ Pqj, simplicity, we first consider a wave packet with A ^ 1 and 6^1. The first 
condition means that it is rapidly oscillating (in the units of the surface gravity) whereas 
the second means that its frequency content is well peaked around the mean frequency A. 
Using the near horizon behavior of Z\, eq. (|5.38 ), the scattered wave packet is 



«$^(uO» = r rfAe^"»^e-^'(^-^)'/V:^ZA(w;/a.ff) (5.42) 

^ ^;^giA(«o+lna„ff+ln(A)/2) g-fe2(A-A)V2 _^ Jf^givr/4^\ 



TT "'-00 L cr, 
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In the second line we have dropped all irrelevant phases and we have used the asymp- 
totic behavior of the F function and the coth, in anticipation to the fact that the main 
contribution of the integral will arise from values of A centered around A. 

By performing the integral by a saddle point approximation (this is perfectly valid 
when 6^ ^ 1), one obtains the following behavior 

exp (- ''" + "";-- + ^'"^'' ) {l + ^/'"fy] . (5.43) 

Hence once Uq> — Ina^g, the image on f = vanishes like a Gaussian in uq. 

This strong dissipation is only valid for tight wave packets in A, i.e. for 6^1. Instead, 
in the opposite regime 6 <^ 1, for tight wave packets in position space, the dissipation is 
milder. Indeed, in the limit 6 — > 0, the Gaussian factor in eq. ( |5.42| ) can be ignored in the 



large Uq limit. Then, the main contribution comes from the first pole of the F function in 
the positive imaginary A axis. In this regime the decrease of the wave packet is given by 

g-(«0+ln o"cff)_ 

In brief, as long as the mean position in i;; on f = of the wave packet is larger 
than AcTeff, its image is unaffected by the metric fluctuations since Zx still behaves as ty*'^. 
Instead, once it enters the near horizon region, its amplitude rapidly decreases. 



6 Conclusion 

It is perhaps appropriate to list what we learn from our analysis. 

1. When the relative width of the smeared horizon is small enough, i.e. when brEul^EU - 
(Jeff -C 1, metric fluctuations in the near horizon geometry affect the asymptotic 



properties of Hawking radiation only slightly, in the second order of a^^, see ([4.17|) 



The reason for this stability can be seen from the short distance behaviour of the 
m-Green function, see ( |4.7|) . Indeed, for |A?/| ^ 1, one recuperates the usual 
Hadamard behavior which guarantees point 1. 

Backward propagated wave packets representing Hawking quanta of energy A are 
dissipated when their Doppler shift frequency id^ reaches cr^~^, i.e. when their sep- 
aration in r from the event horizon approaches Aa^ff, see (p. 421). 



The attentive reader will notice that points 2 and 3 are rather difficult to conciliate. Indeed 
if the behavior of the in-Green function is not modified why for is point 3 relevant? In 
other terms what is the high energy behavior of the theory? Unmodified as suggested by 
point 2 or dramatically modified as indicated by point 3? 

To answer these questions one should first reconsider the domain of validity of the 
scheme we used. Our scheme is based on eq. ( ^.21 ) which represents the free propagation 
of a test-field in a fiuctuating geometry. This is an approximative description which 
completely neglects the back-reaction effects induced by the field itself. In other words 
we are working in the regime when the metric fluctuations are not signiflcantly affected 
by the energy density carried by Therefore, our main hypothesis is that there is an 
intermediate regime in which the metric fluctuations induced by all the other degrees of 
freedom cannot be neglected whereas the back-reaction effects due to 99 can be neglected. 
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For frequencies below this regime, we have proven that the metric fluctuations play no 
role, see eq. (|3.41|) . For higher frequencies, we cannot say much. Nevertheless it is most 
probable that the stochastic description we used fails. Therefore we can trust the behavior 
of the in-Green function only for separations within the intermediate domain. 

Having clarified this point, we can now explain why the high frequency behavior of 
the in-Green function differs so much from that of backward scattered waves. The reason 
is the following. Being a function of the difference in V^, the in-Green function is hardly 
sensitive to the metric fluctuations in the coincidence point limit, see eq. (|4.7| ). On 
the contrary, as emphasized after eq. ( |5.31D , backscattered wave functions defined on 
are sensitive to the metric fluctuations they have encountered when evaluated near 
the horizon. Moreover, since their frequency is blue shifted, they are inevitably strongly 
affected by the metric fluctuations. 

If this explains that the behaviour of the in-Green is perfectly compatible with that 
of backscattered waves, it does not tell us what happens to these waves and their energy 
density when their amplitudes diminish. To describe the fate of their energy density 
one must consider the dynamics of the degrees of freedom which engender the metric 
fluctuations. Indeed, one must go beyond the stochastic treatment of these fluctuations 
in order to be able to describe the trans-Planckian momentum recoils which shall inevitably 
be induced by Hawking photons when they are traced backward near the event horizon, 
T\ IHH] for preliminary attemps to describe this physics. 



see 



In brief, the main outcome of the paper is to have provided physical foundations in 
terms of metric fluctuations to the concept of effective propagation of light near a black 
hole horizon. 

This allows to address in a rational scheme the question of the domain of validity of 
this effective propagation. It also provides an explanation for the vexing question of the 
apparent violation of local Lorentz invariance]^, 0, [1^. The neatest way to characterize 
this violation is to focus on the near horizon behavior of a monochromatic mode 0^"*. In 
the absence of modification of the dispersion relation, this mode behaves as w*^ where 
w = 2r — 1. Hence there is no length which allows one to distinguish low from high 
momenta. This absence is a consequence of the local Lorentz invariance of theories based 
on the usual Dalembertian. On the contrary, when dealing with a modified dispersion 
relation, one breaks this invariance since the new dynamical equation is written in a 
preferred frame. For acoustic black holes this makes good sense since both the frame and 
the critical lenght, which characterizes what "high" frequency means, are given by the 
constituents of the fluid. On the contrary, it is rather unclear to see the origin of such a 
preferred frame for a gravitational black hole. One of the main virtues of the present work 
is to provides a simple answer to this puzzle. Indeed the ensemble of metric fluctuations 
unambiguously determines, a, the constant spread in r (measured along f = 0) of the 
distribution of the backward propagated rays representing the event horizon. Because of 
the hypothesis of stationarity metric fluctuations, the modified equation governing light 
propagation has a simple and stationary expression in the v, r coordinate system. In 
particular, the cut-off lenght a appears only through powers of <jdr\v In this case what 
might be interpreted as the origin of a "violation of Lorentz invariance" originates from 
the ensemble of stationary metric fluctuations. 
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A The leading non-linear corrections 



In this Appendix we shall demonstrate the validity of using the linearized and simplified 
expression, eq. ( p.20|) , of V"^(m) for obtaining the dominant non-linear effects in wo on the 
propagation of the field described by eq. (^.91). 

To simplify the analysis, we consider the Fourier transform with respect to w of the 
image of the pulse defined on J'~ evaluated in a given realization: 



[u\w 



Using eq. 



the mean image is given by 

1 r2n 



A+(m|p) 



(27r)3/2 



the integral over w being trivial. 

The analysis of the null ray propagation in a fluctuating geometry shows 
complete expression for V^{u) is of the form 



V^{u) = -1 - e"" + Fo + Fie"" + 0{e 



-2u\ 



(A.l) 



(A.2) 



that the 



(A.3) 



where the correction terms Fq and Fi are series in wq. The corrections terms multiplied 
by e~" or higher powers of it are not important in the late time regime since they become 
arbitrarily smaller than those of Fq. Hence we can drop Fi. The remaining correction 
term Fq fixes the precise value of the image on J'~ of the exact horizon in every realization 
of the fluctuating geometry. It also determines the late times corrections to A~^{u\p). It 
has the following expansion in powers of wq 



Fq = Wo sin((/) + 0o) + wlfiiu) sin(20 + 0i) + 0{wl) 



(A.4) 



Let us substitute ( |A.4| ) into ( [A. 21 ) and show that the leading non-linear corrections to 
all come from the linear term in wq. 

In the low p regime, i.e. for p ^ I/wq, all correction terms are negligeable since 
Wo ^ 1 in our dimensionless units. Upon reaching the high p regime, for p ~ I/wq, all 
terms of the form pw^ for n > 2 are still very small. Using this fact we can rewrite (|A.2| ) 
as 



A+(m|p) 



g-jp(l+e-") 

" (27r)3/2 

g-ip(l+e-") 

(27r)3/2 ^ 



2n 



^ipwo sin (f> 



1 + ipwlfi{uj) sin(20 + + 
Gq{pwq) + pwIGi{pwq) + pwlGiipwo) + ... . 



(A.5) 



28 



The form of the integrals imphes that functions Gnipwo) are regular. For pwo = they 
are all finite, while for pwo —>■ oo, Gnipwo) are of the same order of magnitude. Therefore 
the corrections to the leading term Gq remain small as long as pw"^ <^ 1. If is of the 
order of the Planck length, the regime of validity of the approximation, p <^ Wq'^, goes 
far beyond the Planckian scale. Having established this result, we can now make more 
precise the neglection of the Fi term in eq. ( [A.3| ) in the late time regime. It is sufficient to 
have e~" < wo, since the linear term in wo in Fi behaves in this regime like the quadratic 
term of Fq. 

In other words, at late times and for all the regimes starting from small p till momenta 
much higher than the Planck one, Go{pwo) (= 27iJo{pwo), see (|3.3CI| )) gives the leading 
non-linear corrections to A~^{u\p). This proves the claim since Gq{pwq) entirely comes 
from the linear term in Fq. 



B Backward Scattering 

In this Appendix we discuss the properties of the backscattered wave of out-frequency A 
when the average is performed over the phase only. 

Its image on J'~ or on f = 0, ^'(a), can be calculated by making a Fourier transform 
as it was done for <^ ^ in Section 5. Here, we shall offer an alternative way to 

compute this image. We start by regularizing the out wave ^°\^{y), see eq. ( ^.24|) . One 
can view this distribution as the limit e — on the real w = y axis of the difference of 
two analytic functions. Explicitly one has 

VA^ifiT'M = i^(w)w'^ = fe"^/2^e - iw)'^ - e-^^/^^e + iwY^) (B.l) 

2 smh(7rA) ^ ^ 

Notice that the out-mode has been expressed by a Bogoliubov relation as the difference of 
two in-modes which are perfectly well defined on the horizon. We can then safely apply 
D to each of them. We shall make use of the relation 



where P^^ ^ denotes the generalized Legendre function. Eq. ( p.2| ) was derived by Filon 
in 1903 and cited in the book by Watson PB[, p. 51, who proposed to the reader to 



prove this formula as an exercise! Using this equation, we obtain 

1 1 



A/47rA 2sinh(7rA) 



[e - iw)"^ + vol 




IW 



iwV + Wn 



where P^{z) is the Legendre function. We have introduced e positive and infinitesimal in 



order to specify the phase of y (e — iwY + in the three sectors: \w\ < Wq and \w\ > wq. 

Before presenting the detailed behavior of $^-) a few simple deductions can be made. 
First, in the limit — > 0, reduces to the usual out wave function given in eq. (|5.24| ) 
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since Pix{l) = 1. Secondly, one finds, as one should do, that vanishes for w < Wq 

since vanished for negative y. Thirdly, the spread of the effective horizon w = ±Wq 

clearly appears in 

Figure ^ illustrates the behavior of for A = 5. This behavior is typical for values 
A greater than 1. For smaller frequencies, the amplitude of the oscillations between ±Wq 
is much smaller than 1. We have presented only the real part of ^^-j since the imaginary 
part behaves similarly. 

The main features of directly follow from its properties near the singular points 
w = oo and w = iwo- The first manifestation of the metric fluctuations shows up around 
w = 20wo where the amplitude of the oscillations start to decrease, see (fig. 5a). This 
modulation of the amplitude arises from the correction term of the asymptotic behavior 
of ^^xy Using eq. 8.1.2 in one has 



The second effect is located for values of w slightly bigger than wq, see (fig. 5b). 
One has a rapidly oscillating function with a decreasing amplitude. The origin of this 
structure can be understood from eq. ( ^.91) : the logarithmicly divergent behavior of ip'^^ 
at w = has been shifted with some coherence to w = Wq because of the stationarity 
of the shift with respect to at its maximal value. This behavior can be seen from the 



analytical property of the function. Using eq. 8.1.5 in [0, one finds that is a sum 
of two hypergeometric functions. The first one controls the overall shape of the function 
near w = Wq whereas the second one provides the rapidly oscillating behaviour. In the 
limit w Wq one has 

where Ai = T^TT-^/^T{iX + l/2)T-\iX + 1) and As = 2-'^-^T{-iX- l/2)T-\iX). Notice 
that one recovers arbitrarily high frequencies in this behavior. Therefore one can fear that 
one has simply shifted the trans-Planckian problem from w = Otow = wq. This is not the 



case for two reasons. First the amplitude of these oscillations decreases as y — Wq for 
w ^ Wq. Therefore the amplitude to find high frequencies decreases. Secondly, when one 
averages over the amplitude Wq, this structure is erased, as shown in Figure 4. Therefore, 
it is an artifact of our simplified averaging procedure. 

The third salient fact is that on the other side of Wq no rapid oscillations are found, 
see (figs. 5b-c). There are nevertheless oscillations but, as indicated in the next equation 
(obtained again from eq. 8.1.5 in |3^) they are suppressed by sinh(7rA) with A = 5. 




*rA)H = ^^^ + TT^I — -1) +0{ — -l]} (B.6) 

The fourth point is that is regular around w = 0, see (fig. 5d). Analytically one 



has 




(B.7) 
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where ii = 7r-i/2r(l + ^)r-i(l + f ) and = - Yvr-i/2r(|i)r-i(i + f ). Therefore, the 
trans-Planckian reservoir of oscillations which was present in <^^"* has been eliminated. 
This is the main physical result. 

The fifth point concerns the structure when w approaches —wq from above. Analjd;i- 
cally, one has 

^7x)M = -T^i-^'r' 1-^^^ ( -2 - iV^'^^' + oK-l]\ (B.8) 

[ tanh(7rA) {w^ J {w^ ' ' 

As in eq. ( [B.5| ), the amplitude of the fluctuations deceases like a square root when 
approaching the boundary w = —Wq. As already mentioned, for smaller values of w one 
finds identically zero. One can apply to this fifth point the remarks made in the second 
point. 

In brief, the main physical result is the fourth point. Moreover, in contradistinction 
with the second, third and fifth points, it is stable if one considers an ensemble of ampli- 
tudes Wo- The physical consequence of the disparition of the reservoir of oscillations for 
w < Wo is that the backward propagation of any wave packet will stop around wo- See 
Section 5.2 for more details. 
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